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ABSTRACT 

This paper addresses available constraints on mass models fitted to rotation curves. Mass models of disk 
galaxies have well-known degeneracies, that prevent a unique mass decomposition. The most notable is due to 
the unknown value of the stellar mass-to-light ratio (the disk-halo degeneracy); even with this known, degen- 
eracies between the halo parameters themselves may prevent an unambiguous determination of the shape of the 
dark halo profile, which includes the inner density slope of the dark matter halo. The latter is often referred to 
as the "cusp-core degeneracy." We explore constraints on the disk and halo parameters and apply these to four 
mock and six observed disk galaxies with high resolution and extended rotation curves. Our full set of con- 
straints consists of mass-to-light (M/L) ratios from stellar population synthesis models based on B — R colors, 
constraints on halo parameters from A/-body simulations, and constraining the halo virial velocity to be less 
than the maximum observed velocity. These constraints are only partially successful in lifting the cusp-core 
degeneracy. The effect of adiabatic contraction of the halo by the disk is to steepen cores into cusps and reduce 
the best-fit halo concentration and M/L values (often significantly). We also discuss the effect of disk thick- 
ness, halo flattening, distance errors, and rotation curve error values on mass modeling. Increasing the imposed 
minimum rotation curve error from typically low, underestimated values to more realistic estimates decreases 
the x 2 substantially and makes distinguishing between a cuspy or cored halo profile even more difficult. In 
spite of the degeneracies and uncertainties present, our constrained mass modeling favors sub-maximal disks 
(i.e., a dominant halo) at 2.2 disk scale lengths, with Vdi S k/Vtot ^ 0.6. This result holds for both the un-barred 
and weakly barred galaxies in our sample. 

Subject headings: dark matter — galaxies: fundamental parameters — galaxies: halos — galaxies: kinematics 
and dynamics — galaxies: spiral — galaxies: structure 



1. INTRODUCTION 

There has been significant debate recently about the shape 
of dark matter density profiles, especially regarding their in- 
ner slope, a. 1 Based on cosmological A/-body simulations 
(Navarro, Frenk, & White 1996; Navarro, Frenk, & White 
1997; hereafter NFW), the dark matter halo profile appears to 
be independent of mass and has an inner logarithmic slope 
a = 1. More recent, higher resolution simulations suggest 
that the density profiles do not converge to a single power- 
law at small radii. At the smallest resolved scales (~ 0.5% 
of the virial radius) profiles usually have slopes between 1 
and 1.5 (Moore et al. 1999; Ghigna et al. 2000; Jing & Suto 
2000; Fukushige & Makino 2001; Klypin et al. 2001; Power 
et al. 2003; Navarro et al. 2004; Diemand et al. 2004). 

' Values of a range from (core) to 1.5 (cuspy). We define the dark matter 
profile and a in Equation (4). 



At large radii all simulations find density profiles with 
slopes a ~ -3, which is inconsistent with the isothermal 
(p oc r~ 2 ) profile. The determination of the dark halo slope 
based on mapping the outer density profile of galaxies is dif- 
ficult, owing mainly to a lack of mass tracers at large radii. 
Prada et al. (2003) find that the line-of-sight velocity disper- 
sion of satellite galaxies declines with distance to the primary, 
in agreement with a p oc r" 3 density profile at large radii. 

The determination of a based on data at smaller radii is 
complicated by the unknown value of the stellar mass-to-light 
ratio, Tj. This has led to dedicated analyses on dwarf 2 and 
low surface brightness 3 (LSB) galaxies that are believed to be 

2 Dwarf spiral galaxies are usually defined as having a maximum rotation 
velocity v m ax < 100 km s and/or a total magnitude Mb > —18. 

3 An LSB galaxy is usually defined as a disk galaxy with an extrapolated 
central disk surface brightness roughly 2 magarcsec -2 fainter than the 
typical value for HSB galaxies of p® = 21.65 (Freeman 1970). 
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dark matter dominated at all radii (de Blok & McGaugh 1997; 
Verheijen 1997; Swaters 1999). 

It has been suggested that rotation curves of dwarf and LSB 
galaxies rise less steeply than predicted by numerical simula- 
tions based on the cold dark matter (CDM) paradigm (Moore 
1994; Flores & Primack 1994; de Blok & McGaugh 1997; 
McGaugh & de Blok 1998; de Blok et al. 2001a,b). How- 
ever, a number of observational uncertainties cast doubt over 
these early conflicting claims. These include beam smearing 
for Hi rotation curves (Swaters et al. 2000; van den Bosch 
et al. 2000), high inclination angles and Ha long-slit align- 
ment error (Swaters et al. 2003a), and non-circular motions 
close to the center of galaxies (Swaters et al. 2003b). Many of 
these uncertainties can be quantified or eliminated by measur- 
ing high-resolution two-dimensional velocity fields (Barnes, 
Sellwood, & Kosowsky 2004). At optical wavelengths, these 
can be obtained via Fabry-Perot interferometry (e.g., Blais- 
Ouellette et al. 1999) or integral field spectroscopy (e.g., An- 
dersen & Bershady 2003; Courteau et al. 2003). 

Despite a low ratio of baryonic to non-baryonic matter in 
dwarf and LSB galaxies, practical limitations in accurately 
determining the circular velocity profile have prevented a reli- 
able determination of the dark matter density profile for those 
galaxies. Furthermore, the predictions of numerical simu- 
lations are weakest on the (small) scales of dwarf and LSB 
galaxies. By comparison, for high surface brightness (HSB) 
galaxies the kinematics is easier to measure and the expected 
dark halos can be better resolved in numerical simulations, but 
the more prominent stellar component often hinders a unique 
mass decomposition. 

In principle, if the disk mass-to-light ratio, T<j, and the 
gaseous mass distribution are known, the contribution from 
the dark halo to the overall potential can be determined. How- 
ever, extracting the parametrized halo profile with this proce- 
dure is complicated owing to a degeneracy between the halo 
parameters themselves (e.g., van den Bosch & Swaters 2001). 
Furthermore, various evolutionary processes may alter the 
dark halo density profile from that found in dark matter-only 
simulations. The dissipation of the disk is thought to com- 
press the dark halo distribution through adiabatic contraction 
(Blumenthal et al. 1986; Flores et al. 1993), while other pro- 
cesses such as feedback, mergers, spin segregation (Mailer & 
Dekel 2002; Dekel et al. 2003), pre-processing of dark halos 
(Mo & Mao 2003), and bar-driven dark halo evolution (Wein- 
berg & Katz 2002) are thought to lower the concentration and 
central cusp of dark matter halos. 

In this paper we discuss and apply mass modeling con- 
straints in an attempt to break internal modeling degeneracies 
and thus determine the best parameterization of the dark halo. 
We present our mass models in §2 and their degeneracies in 
§3. The mass model constraints are presented in §4. We then 
apply these constraints to six galaxies from Blais-Ouellette 
(2000). The data are presented in §5, and the models are ap- 
plied to the data in §6. The effects of rotation curve errors, 
distance, disk thickness, and halo flattening are discussed in 
§7, and a summary is offered in §8. 

Throughout this paper r and R refer to the radius from the 
galaxy center in spherical and cylindrical coordinates, respec- 
tively. Whenever necessary, we also adopt a value of the Hub- 
ble constant 4 H given by h = H /100 = 0.7. 



The current best estimate of the Hubble constant is H 
8kms-' Mpc~' (HST H Key Project; Freedman et al. 2001). 
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2. MASS MODELS 

Our mass models consist of three main components for each 
disk galaxy: a thick stellar disk (hereafter the "disk"), an in- 
finitesimally thin gas disk (hereafter the "gas"), and an oblate 
dark halo (hereafter the "halo"). In general, disk galaxies may 
also have a bulge component, although for simplicity we limit 
our analysis to nearly bulge-less systems. Assuming that the 
matter distribution is axially symmetric and in virial equilib- 
rium, the total circular velocity is given by 



(1) 



at each radius R. Each of the three components is described in 
more detail below. We compute the circular velocities of the 
disk and gas using formula A. 17 of Casertano (1983). The 
best-fitting mass model is determined by fitting V c j rc to the 
observed circular rotation velocity, V rot , by minimizing the \ 2 
-statistic with a non-linear optimization scheme. 

2.1. Stellar Disk 

We model the disk with the following density profile (van 
der Kruit & Searle 1981): 



£(fl)sech 2 (z/z ) 
2zo 



(2) 



where £(/?) is the disk surface density profile and zq is the 
vertical scale height. Unless otherwise stated, we compute 
£(/?) from the observed surface brightness profile. 

The vertical scale height is parameterized in terms of the in- 
trinsic disk thickness, = zo/Rd, where /?d is the disk scale 
length. Unless otherwise stated, we adopt q^ = 0.25 (Kregel 
et al. 2002; Bizyaev & Mitronova 2002). We explore the ef- 
fect of disk thickness in §7. 

2.2. Gaseous Disk 
We model the gas disk with the following density profile: 

Ps - dS (R,z) = 6(z)Z m (R)/fm (3) 

where S(z) is the Kronecker delta function, £m(/?) is the sur- 
face density of neutral hydrogen, and /hi is the fraction of gas 
in Hi. We adopt f m = 0.75 (e.g., Blais-Ouellette et al. 2001); 
other authors take 0.71 < / ffl < 0.77, although the exact value 
is not critical. 

Some spiral galaxies show a central depression in the Hi 
density, likely due to the gas being present in a different form 
(ionized or molecular) and/or partially or completely con- 
sumed in previous episodes of star formation. A central de- 
pression and hence a positive radial density gradient result in 
an outward radial force or negative V g 2 as . We represent this 
as negative velocity on the gaseous component of the rotation 
curve. 

2.3. Dark Halo 

To take account of the uncertainties in the predicted inner 
halo density profiles and to allow for halos with flat central 
density profiles while preserving the p oc r" 3 dependence at 
large radii, we use the following density profile (hereafter 
ALP profile; Kravtsov et al. 1998): 



PhaloM ; 



P0 



(r/r s r(l + r/r s ) 



3-a " 



(4) 



This density profile has an inner logarithmic slope of -a and 
an outer logarithmic slope of -3. For a = 1 this reduces to 
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log(r/r_ 2 ) 

FIG. 1 . — Logarithmic slope of the 'ALP' halo density profiles for a = 0, 1, 
and 1.5 andc= 10 plotted from 0.4% rioo < r < T2oo- The highest resolution 
N-body simulations can resolve the density profile over this range. Typi- 
cal density profiles lie between the a = 1 and a = 1.5 lines. For compari- 
son with the 'ALP' parameterization we show the fitting function of Navarro 
et al. 2004. 
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FIG. 3. — Circular velocity profiles for halos in Fig. 1 normalized to the 
virial velocity, V200 The bottom panel shows the differences with respect 
to the o = l profile. Note that the differences between the profiles are most 
conspicuous for radii less than ~ 0.5 r_2- Also shown is a halo with a = 1.2 
(red dot-dashed line), which is effectively indistinguishable from the Navarro 
et al. (2004) fitting function. 



the NFW profile, and at the scale radius, r s , the slope of the 
density profile is -2. 

However, for different values of a, r s corresponds to differ- 
ent density slopes. To enable an easier comparison of scale 
radii, we replace r s with r_2, the radius where the slope of the 
density profile is -2. With the conversion r_2 = (2- a)r s . 

Figures 1, 2, and 3 show the logarithmic density slopes, 
density profiles, and circular velocity profiles for halos with 
a = 0,1, and 1 .5 and the fitting function from Navarro 
et al. 2004 (which has effectively a ~ 1.2). 
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2.3.1. Oblate/Prolate Density Profiles 

Typically rotation curve analyses assume a spherical dark 
halo even though CDM simulations suggest triaxial shapes 
for collapsed structures, with typical axis ratios c/a = 0.5-0.7 
and b/a = 0.7-0.9 (Dubinski & Carlberg 1991; Jing & Suto 
2002; Tinker & Ryden 2002). However, the dissipative infall 
of gas in non-baryonic dark halos suppresses triaxial struc- 
tures, leading to halos with an oblate shape (Katz & Gunn 
1 99 1 ; Dubinski 1 994; although further investigation is needed 
to quantify this effect). 

This tentative conclusion agrees with a variety of observa- 
tions that find axially symmetric disks, with eccentricity e < 
0.05 (b/a > 0.9987; Combes 2002 and references therein). 
The flattening of the halo is not easily measured, but vari- 
ous techniques, including the flaring of Hi disks, polar rings 
around spiral galaxies, and X-ray isophotes of elliptical galax- 
ies, suggest oblate halos with an axis ratio, q = c/a, ranging 
from 0.1 to 0.9 (Combes 2002). 

Thus, we are compelled to study the effects of axially sym- 
metric dark halos (b/a = 1) in our mass models. We gener- 
alize the density profiles to the family of axially symmetric 
ellipsoids by setting p(r) = p(m), where m 2 = R 2 + z 2 /q 2 - 

2.3.2. c — V200 parameterization 

We choose to parameterize the density profile by the cir- 
cular velocity at the virial radius, V wu , and the concentration 
parameter, c-2 = R v i T lr-2- Here R wu is the virial radius in the 
z = plane. We choose to define the virial radius, R wu , as the 
radius where the mean density of the halo is A v ; r times the 
critical density, 



FIG. 2. — Density profiles for halos in Fig. 1 normalized to the density at 
the scale radius, p_2 ■ The bottom panel shows the differences with respect to 
the a = 1 profile. 



M(m vir ) 

P(fWvir) = 4 5- = A v 



■Pcrit, 



(5) 
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FIG. 4. — Effect of halo flattening q on circular velocity profiles for a = 1 
halos with c = 10. The bottom panel shows the fractional differences with 
respect to the q = 1 case. Note that these differences increase with decreasing 
radii. 



where p c „ t = is the critical density of the universe. With 
these definitions m v i r , and hence R v [ r , will be invariant under 
changes of q. Unless otherwise stated, we adopt A v ; r = 200, 
although in currently favored ACDM cosmologies, at redshift 
zero, the virial radius occurs at A™ = 337 SIm — 100 (Bryan & 
Norman 1998). The choice of A v i r does not affect the density 
profile, but the virial radius changes by a factor of <~ 1.3. 

With these definitions and using Equation 2.91 in Binney 
& Tremaine (1987) for the computation of V c { rc , the velocity 
contribution of the halo specified by V200. c -2, ce, and q is 
given by 



y^ lo (x,z = 0) = y 2 2 00 



H(x,a,q)/x 
fi(c- 2 ,a,q)/c-2 



x = R/r. 2 (6) 



where 



fji(x,a,q) = 



a-3 



y 2 -"[l-(2-q)y] 
y/\-(\-q 2 )f/x 2 



dy. 



(7) 



With the above definitions we can express the relationship 
between R vn and V v ; r as 



Kir \ 2 _ h 2 ( Avir \ fJ,(C-2,a,q) 

R n J ~ \200j m(c-2,«,1) 



(8) 



with V v ; r and R v i r in kms 1 and kpc, respectively, and h = 
#o/100 = 0.7. 

Fig. 4 shows the effect of q on the circular velocity of the 
halo, normalized by V v ; r . Note that for a given R v \ r , c_2, and 
a, y v ; r increases as q decreases. Oblate (q < 1) halos result 
in higher circular velocities, especially near the center, while 
prolate (q > 1) halos result in lower circular velocities. 

2.3.3. Adiabatic Contraction 

We model the response of the dissipation-less halo to the 
infall of the dissipational baryons as they cool and settle 
into a disk following Blumenthal et al. (1986) and Flores 
et al. (1993). This assumes that the collapse of the baryons 



is slow, the matter distribution is spherically symmetric, and 
particles move on circular orbits. Then the adiabatic invariant 
is simply rM(r), where M(r) is the mass enclosed by radius 
r. With the further assumptions that the dark matter parti- 
cles do not cross orbits, M DM (?"f) = M DM (nX where r\ and rf 
are the initial and final radii of the disk, respectively, and that 
the baryons are initially mixed with the dark matter with a 
baryon fraction fg = Mb/ (Mb + Mom), then given the initial 
dark halo distribution Mdm^;) and final baryonic mass distri- 
bution Mfi(rf), the final radius, rf, can be obtained by solving 

(9) 



rf[M B (r f )+M DM (n)] = r 1 M DM (r i )/(l -/ B ). 



We obtain the baryonic mass from the observations of stars 
and gas in the disk (§2.1), for a given Td and distance, and 
assume that the fraction of baryons in the halo is negligible. 
For the mass of the halo we assume the virial mass, M v ; r . 

The effect of adiabatic contraction on the density and cir- 
cular velocity distributions can be quite substantial. We illus- 
trate this effect in Figure 5. This shows the circular velocity 
of the initial and final halo and final disk. Here the disk is 
exponential with Rd = 2kpc, p,* = 20 mag arcsec" 2 (R-band), 
and T* = 1.0 and 0.25. The effect of adiabatic contraction is 
largest for halos with low values of a and c_2, such that halos 
with initial cores end up with cusps. For very cuspy halos, 
the halo can expand in the very center, as the final baryonic 
mass within r; is less than the initial baryonic mass within r\. 
Note that although the mass is more centrally concentrated 
after adiabatic contraction, the formal concentration, c_2, can 
stay the same. 

Despite the simplifying assumptions, the validity of the adi- 
abatic approximation of Blumenthal et al. (1986) has been 
confirmed down to 10~ 2 r_2 in a study of the response of a 
dark matter halo to the growth of an exponential disk in high- 
resolution .A-body simulations (Jesseit et al. 2002). However, 
Wilson (2003) and Gnedin et al. (2004) claim that under more 
general conditions the standard model for adiabatic contrac- 
tion systematically overpredicts the contraction in the inner- 
most regions, while slightly underpredicting the contraction 
at larger radii. Therefore, we use the standard model for adi- 
abatic contraction to provide an upper limit on the effect of 
adiabatic contraction. 

3. MODEL DEGENERACIES 

Several degeneracies exist between the model parameters, 
which may prevent a unique mass decomposition. These can 
be divided into "disk-halo" and "cusp-core" degeneracies. 

The "disk-halo" degeneracy occurs between Td and the 
halo parameters. Equally good fits, in a x 2 -statistical sense, 
where x\ is the reduced-% 2 , can be obtained with a wide 
range of Td from zero to a maximum disk (e.g., van Albada 
et al. 1985). To break this degeneracy, we need a priori knowl- 
edge of Td, or Vdisk/Vtot- If the disk thickness and halo flatten- 
ing are ignored, this fixes the density profile of the dark halo. 
However, in practice, with errors on the circular velocities of a 
few kilometers per second finite spatial resolution, and limited 
extent of the rotation curve, degeneracies between the halo pa- 
rameters themselves often prevent a unique parameterization 
of the halo density profile (van den Bosch & Swaters 2001). 
To break this degeneracy, constraints need to be placed on c_2 
and V200 as well. 

We illustrate these degeneracies with mock rotation curves. 
Our mock galaxies consist of an exponential disk specified by 
= 20 mag arcsec" 2 , Rd = 2 kpc, and =1 and an adiabat- 
ically contracted dark halo with initial parameters c_2 =10, 
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FIG. 5. — Effect of adiabatic contraction on circular velocity and halo density slopes for halos with a =0, 1.0, and 1.5 (left to right), c_2 = 5 (top) and 10 
(bottom), for exponential disks with Ri = 2 kpc and So = 258 and 52 Mq pc~ 2 . For each disk-halo system we show the initial halo (red dotted line), final halo 
after adiabatic contraction (AC; red sold line), disk (blue dashed), and final total circular velocity (black thick solid). 



V200 = 100 kms" 1 , and a = 0, 0.5, 1, and 1.5. We then sample 
the rotation curve in 3" bins up to 2Rd and in 15" bins up to 
8/?^ to simulate Ha and Hi data respectively. We then add a 
random Gaussian error (with a = 4kms _1 ) and assign a con- 
servative measurement error of 4 kms -1 to each data point. 

We fit for c-2 and V200 on a grid of ol and T*, with and with- 
out adiabatic contraction. The results of these fits are shown 
in Fig. 6. The disk-halo and cusp-core degeneracies exist for 
all input values of a and are strongest for a = 1 halos. Thus, 
for these model galaxies, without constraints it is impossible 
to determine or a based on the value alone. When 
fitting without adiabatic contraction, a wider range of T* is 



permitted, including maximum disks for a = halos. We also 
see that the form of the c-2 - a relation is the same for all fits, 
but the normalization is lower for fits with adiabatic contrac- 
tion and a higher input a. 

In order to achieve reliable results out of the mass modeling 
exercise, we must therefore consider independent constraints, 
which we discuss below. 

4. CONSTRAINTS 

4. 1 . Stellar Population Synthesis Models 

Stellar population synthesis (SPS) models can be used to 
place constraints on Tj. The combination of optical and in- 
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FIG. 6. — Best-fitting halo parameters and Xt for mock rotation curves vs. a. All input models have adiabatically contracted (AC) halos with c_2 = 10, 
V200 = 100, and exponential disks with = 20 and = 1.0. The only difference is in the central density slope, a. The input models are indicated by the 
filled circles and have xl = 1 -0. The top panels show fits with AC, while the bottom panels show fits without. The different lines correspond to different fitted 
mass-to-light ratios as indicated. The horizontal line indicates the maximum "observed" circular velocity. Note that we have imposed c_2 > 2. 
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frared photometry, with SPS models, yields T d values accu- 
rate to ~ 40% (Bell & de Jong 2001). The slope of T d versus 
color is fairly independent of the initial mass function (IMF) 
and star formation (SF) history. That slope is also smaller in 
the TT-band than in the B-band, although the zero point of the 
color- Yd relation is itself very sensitive to the IMF. The cal- 
ibration of the Td -color relation by Bell & de Jong (2001) 
relies on the assumption of some galaxies being close to max- 
imal disks, and their values of T d are thus upper limits. 

4.2. Evidence for Sub-Maximal Disks 

A conventional hypothesis to determine an upper limit to 
Td is that disks should be maximal 5 (Carignan & Freeman 
1985; van Albada & Sancisi 1986). This approach works 
well in practice for most HSB galaxies, but dark matter is still 
needed to explain the outer rotation curves of HSB and LSB 
galaxies at almost all radii. The fact that maximum disks can 
match inner rotation curves of HSB galaxies is more telling 
about the degeneracies in mass modeling than the validity of 
the hypothesis itself (Broeils & Courteau 1997; Courteau & 
Rix 1999). Furthermore, Ha rotation curves alone can often 
be fitted by pure disk or pure halo models and thus lack any 
constraining power without the addition of an extended Hi ro- 
tation curve (Buchhorn 1992; Broeils & Courteau 1997). 

By contrast to the maximal disk hypothesis, a variety of 
methods, which are described below, suggest that on average 
HSB disks are sub-maximal with 

(WV t ot) 2 . 2 ^0.6. (10) 

Note that a galaxy with a sub-maximal disk at 2.2R^ can still 
be baryon dominated at 2.2R^ if there is a significant bulge 
component. 

4.2.1. TFR Residuals 

Courteau & Rix (1999) have suggested that sub-maximal 
disks should be invoked to explain the surface brightness in- 
dependence of the Tully-Fisher relation (TFR); they find that, 
on average, HSB galaxies have (Vdisk/Vobs) 2 2 ^ 0.6 ±0.1 
(see also Courteau et al. 2005). Their argument only depends 
on the assumptions that the scatter in the TFR and the size- 
luminosity relation (SLR) is dominated by a dependence in 
Rd and that dark halos respond adiabatically to the formation 
of the disk. 

4.2.2. Velocity Dispersion Measurements 

The peak circular velocity of an isolated exponential disk 
can be related to the vertical velocity dispersion 6 and the in- 
trinsic thickness of the disk, zo, via (Bottema 1993) 

V^ = 0M(Vf) l ,l'ff. (11) 

V ~0 

The factor 0.88 applies to a disk of zero thickness; for a 
thicker disk the peak velocity will be lower. 

In practice, these measurements are difficult since the scale 
height, zo, and scale length, Rd, of the disk cannot be mea- 
sured simultaneously and the vertical velocity dispersions are 
easiest to measure in face-on galaxies although the disk kine- 
matics is hard to determine. The prospects for this method 

5 We adopt the definition of a maximal disk as one that supplies 85 ± 10% 
of the total velocity at 2.2R d (Sackett 1997). 

6 The correction to the velocity dispersion of a disk embedded in a dark 
matter halo is usually negligible (Bottema 1993). 



are improving with the ability to reliably determine the incli- 
nations for low-inclination galaxies using integral field spec- 
troscopy (Verheijen et al. 2004). 

Bottema (1993) found, for stellar kinematic measurements 
in 12 spiral galaxies (with V ™ x > lOOkms -1 ), that more mas- 
sive spirals have larger velocity dispersions with the corre- 
lation (v2)J/ = 2 = {v£)^ Rd = (0.30 ± 0.06) V™?. Substituting 
this into equation (11) and taking the intrinsic disk scale ratio 
Rd/zo = 4.2 ± 1 .5 (Kregel et al. 2002) yields 

V d ™k X /V o r = 0.5±0.2. (12) 

By comparison, Bottema (1993) obtained a mean value of 
0.63, using Rd/zo = 6. 

4.2.3. Gravitational Lensing 

In some rare cases in which a quasar is lensed by a fore- 
ground galaxy and gravitational lensing can be used to place 
an extra constraint on the mass profile, the dynamical analysis 
strongly favors sub-maximal disks (Mailer et al. 2000; Trott 
& Webster 2002). 

4.2.4. Bars and Spiral Structure 

It is generally thought that dynamical friction between the 
bar and halo will slow down the pattern speed of the bar; 
thus, fast bars imply maximal disks (Weinberg 1985; Hern- 
quist & Weinberg 1992; Debattista & Sellwood 2000). How- 
ever, other authors claim that the efficiency of bar slow-down 
by dynamical friction has been overestimated (Valenzuela & 
Klypin 2003) and that the bar pattern speed is not a reliable in- 
dicator of disk-to-dark matter ratio (Athanassoula2003). Cur- 
rent observational data favor fast bar pattern speeds; however, 
only a handful of galaxies have reliable measurements (e.g., 
Debattista & Williams 2004), and most observations are of 
SB0 spiral galaxies, whose large bulge components and red 
colors are consistent with being baryon dominated. 

In terms of late-type spiral galaxies, Weiner et al. (2001) 
modeled the strong shocks and non-circular motions in the 
observed gas flow and find that a high T d , corresponding to 
80%-95% of Vdisk, and a fast-rotating bar are highly favored. 
On the other hand, modeling of the spiral arm structure of 
a few grand-design galaxies by Kranz et al. (2003) yields a 
wide range of Vdi s k/Vtot, from closely maximal to 0.6. It should 
be noted that both these methods are model dependent and 
that while the basic dynamics of bars and spiral structure are 
understood, there are issues that remain to be resolved. 

Courteau et al. (2003) showed that barred and non-barred 
galaxies belong to the same TFR. Thus, if the argument by 
CR99 is correct, barred galaxies would, on average, harbor 
sub-maximal disks. Note that this is consistent with the above 
observations, if there is significant scatter in Vdisk/Vtot, or ^ 
Vdisk/^tot increases with surface brightness. 

4.3. Constraints on V 2 oo 

As shown in Figure 6, the fitted V 2 oo often exceeds the max- 
imum rotation velocity of the galaxy; by restricting V 2 oo < 
Vmax, the parameter space is reduced. For observed galaxies 
we expect > V max provided that the rotation curve flat- 
tens out or declines at large radii, as is typical for extended 
rotation curves of spiral galaxies (Casertano & van Gorkom 
1991). 

The combined analysis of galaxy-galaxy lensing from the 
Sloan Digital Sky Survey (SDSS) and the TFR led Seljak 
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TABLE 1 

Galaxy parameters 



Galaxy Mb Band jUg' c Rj D Vq i 1/ cos(i') V mll x Rn a Rm References 

(mag) (magarcsec -2 ) (kpc) (Mpc) (kms -1 ) (deg) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



NGC3109 -16.35 B 22.3 1.3 1.36 c 403 75±5 3.86 67 2.1 5.0 1,2,1,3 

IC2574 -16.77 R 23.0 2.3 4.0 b 57 75 ±7 3.86 67 4.6 4,2,4,5 

UGC2259 -17.03 r 21.7 1.6 10.3 h 583 41 ±3 1.33 90 1.4 4.8 6,7,8,7 

NGC5585 -17.50 R 21.1 1.9 8.7 b 305 53 ± 1 1.66 92 2.7 7.1 9,10,9,11 

NGC2403 -19.50 r 20.4 1.8 3.22 c 131 60±2 2.00 136 2.1 10.8 12,7,8,13 

NGC3198 -19.90 r 21.0 3.7 13. 8 C 663 72±2 3.24 157 1.8 12.0 12,14,8,13 



REFERENCES. — (1): Jobin & Carignan (1990); (2): Blais-Ouellette et al. (2001); (3): Musella et al. (1997); (4): Martimbeau et al. (1994); (5): 
Karachentsev et al. (2002); (6): Carignan et al.(1988); (7): Blais-Ouellette et al. (2004); (8): Kent (1987); (9): Cote et al. (1991). (10): Blais-Ouellette et al. 
(1999); (11): Drozdovsky et al. (2000); (12): Begeman (1987); (13): Freedman et al. (2001); (14): Corradi et al. (1991); 

NOTE. — Col. (1): Galaxy name. Col. (2): Absolute B magnitude (Blais-Ouellette 2000). Col. (3): Photometry band; where necessary we convert to 
Cousins R assuming B—R = 0.9 or i — R = 0.35 (Jorgensen 1994). Col. (4): Central surface brightness in R band from a fit to the surface brightness profile 
with a marked disk, corrected for inclination (Col. [8]) and Galactic extinction (Schlegel et al. 1998). Col. (5): Scale length of the disk, from a fit to the 
i?-band surface brightness profile with a marked disk. Col. (6): Adopted distance, with distance indicators "c" for Cepheid, "b" for brightest stars, and "h" for 
Hubble distance (h = 0.7) corrected for Virgo-centric flow (150 km/s). Col. (7): Systemic velocity; see HI reference. Col. (8): Inclination; see HI reference. 
Col. (9): Conversion factor we use to correct our Tj for inclination, assuming a thin disk with no extinction. Col. (10): Maximum rotation velocity; see HI 
reference. Col. (11): Radius out to which we use the Ha rotation curve (in Ra). Col. (12): Radius out to which we use the HI rotation curve (in Sj). Col. 
(13): References, in the following sequence: HI, Ha, photometry, D. 



(2002) to postulate that the rotation velocity of early- and late- 
type ~ L* galaxies decreases significantly from its peak value 
at the optical radius (3.2/?d) to the virial radius R20Q, with 

V ^/V(R 20 o)~l.8, (13) 

and a 2a lower limit of 1.4. This implies that the rotation 
curve declines at large radii and is thus inconsistent with an 
isothermal profile (p oc r~ 2 ), unless the velocity excess over 
V(/?2oo) is due entirely to the stellar component, which is un- 
likely for late-type galaxies. This result is consistent with 
Prada et al. (2003), who find p oc r" 3 at large radii. 

Combining this constraint with the above evidence for sub- 
maximal disks, we get the interesting result that the total ve- 
locity at the virial radius is approximately equal to the peak 
velocity of the disk, 

VS~V(R 2Q0 )^V 20Q . (14) 

4.4. Halo Concentration Parameter 

Cosmological simulations suggest a correlation between 
c-2 and V v i r , such that more massive halos have lower con- 
centrations (e.g., Bullock et al. 2001; Eke et al. 2001; Wech- 
sler et al. 2002; Zhao et al. 2003). This is because halos 
with smaller masses collapse earlier, when the universe has 
a higher mean density. In the standard A CDM cosmology 
((is = 1 .0, = 0.3) for V200 =100 kms -1 the mean c_2 ~ 12, 
but with a significant scatter of A log c_2 = 0. 14 (Wechsler etal 
2002). The mean c_2 is roughly proportional to erg and is also 
weakly dependent on S! M such that lower £1 M gives lower c_2- 
Given the current observational uncertainties in these param- 
eters, the mean c_2 could easily be lowered by 25%. In the 
model of Bullock et al. (2001) for 40 < V200 < 160, the 2a 
range in concentration is 6 < c_2 ^ 30. 

5. DATA 

We now apply these constraints to a test sample of disk 
galaxies from Blais-Ouellette (2000). This sample is one 
of the few with rotation curves derived from both two- 
dimensional Ha and Hi velocity fields. The Ha data are 
needed to probe the inner rising part of the rotation curve, 



where resolution effects often limit the reliablity of Hi data. 
However, Hi rotation curves are still needed to probe the outer 
part of the rotation curve (typically extending to twice the op- 
tical radius) and are essential to constrain the halo mass. The 
full Blais-Ouellette sample contains 10 galaxies, with a wide 
range of luminosities and surface brightnesses. We restrict 
our analysis to six of these galaxies that do not have signifi- 
cant bulges. Accurately decomposing the disk and bulge com- 
ponents and determining radial Yd gradients requires near-IR 
imaging (e.g. K band) and optical color profiles. We will re- 
turn to these issues in a forthcoming paper. Table 1 gives the 
optical and kinematic parameters of the sample galaxies, as 
well as references for the data sources. While this sample is 
by no means complete, it provides a representative selection 
of LSB and HSB "bulge-less" disk galaxies against which we 
can test our constraints. 

5.1. Rotation curve errors 

Ideally, the rotation curve errors would be normally dis- 
tributed and indicate the uncertainty in the circular velocity 
at a given radius. In practice, the errors are often defined in 
some ad hoc fashion, such as assigning a constant value to 
each velocity bin, and cannot be used to place confidence lev- 
els on fitted model parameters. We assume that the observed 
rotation velocity is equal to the circular velocity and there- 
fore that non-circular motions (e.g., bars, streaming motions, 
pressure support) are not significant. 

Blais-Ouellette (2000) computes errors on the Ha rotation 
curves using a j\f~N in each ring of a tilted ring fit to the ve- 
locity field. By contrast, the errors on the Hi rotation curves 
in our sample are computed using the velocity difference be- 
tween the approaching and receding sides of the galaxy. To 
be consistent, we re-compute the errors on the Ha rotation 
curves by taking the maximum of a/ \/N and the velocity dif- 
ference between the two sides weighted by the number of 
points on each side. We also impose a minimum error of 
2km s -1 (to be consistent with the Hlobservations), although 
the actual uncertainty is probably larger. We investigate the 
effect of a larger minimum error value in §8. 
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FIG. 7. — Best-fitting halo parameters and x? vs. a for a range of (as specified in the figures). We show fits with (AC) and without (NAC) adiabatic 
contraction and with (+) and without constraints on , c_ 2 , and V 2 00- The dashed horizontal line indicates the constraint on V 2 oo, from the maximum observed 
velocity; for NGC 2403 and NGC 3198 we also show Vmax/1.4 and Vmax/1.8. For NGC 3109 and IC 2574 the rotation curves are still rising at the last measured 
point; here we set V200 < 100. 
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FIG. 7. — Continued. 
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FIG. 8. — Sample rotation curve decompositions for all six galaxies. Shown are best fits (minimum \ 2 ) f° r ce = and a = 1, with (AC) and without (NAC) 
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FIG. 9. — Effect on best-fit parameters of adiabatic contraction for halos with a = (red) and a = 1 (black). Dashed lines show fits without AC, while solid 
lines show fits with AC. 
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5.2. Photometry uncertainties 

In principle, observed surface brightness profiles should be 
corrected for projection and extinction effects to bring the 
profiles to their intrinsic face-on values. These have oppo- 
site effects on the surface brightness profile (the former leads 
to an over-estimate of the surface brightness, while the latter 
causes it to be under-estimated). In practice, these correc- 
tions are fraught with significant uncertainties, especially for 
extinction correction. For simplicity we only correct for in- 
clination assuming that the disk is optically thin and of zero 
thickness. Thus, the observed surface density will be a fac- 
tor of 1 / cos(;') greater than the face-on case, where i is the 
inclination angle. The adopted inclination angles are given 
in Table 1 . We correct for Galactic extinction using the red- 
dening values of Schlegel et al. (1998). The effect of internal 
dust extinction is to reduce the observed surface brightness 
from its intrinsic (stellar) value and is expected to have some 
radial dependence. Thus, the simple dust-free model that we 
adopt, all other variables being equal, gives an upper limit to 
the T* 

6. ROTATION CURVE FITS 

As a result of certain covariances between model param- 
eters, we fit for c_2 and V200 on a grid of and a to pre- 
vent our fitting routine from getting trapped in local minima. 
The results of these fits are shown in Figure 7, where we plot 
the reduced x 2 an d best-fit c_2 and V200 against a for a range 
of T* from zero to maximum disks. Examples of rotation 
curve decompositions for a = and a = 1 halos are shown in 
Fig. 8. Clearly, without constraints the models are highly de- 
generate, with acceptable fits being possible covering a wide, 
ill-constrained range of a and T*. 

6. 1. Effect of adiabatic contraction 

In Figure 9 we plot the Xr an d best fit c_2 and V200 versus 
, for a = and a = 1 halos. Adiabatic contraction causes 
the rotation velocity of the halo to rise, especially at small 
radii. Thus, to obtain a comparable circular velocity profile, 
the concentration of the pre-contracted halo must be lower. 
As T J increases, this effect becomes more significant and can 
substantially alter the halo parameters, even for sub-maximal 
disks. For close to maximum disks acceptable fits cannot be 
obtained without breaking the weak constraint: c_2 > 2. Even 
after exclusion of the high Vdi s k /V tot solutions, the disk-halo 
and cusp-core degeneracies still remain. In general, the best- 
fit Xi values are lower with adiabatic contraction than without. 
In cases in which a maximum disk is the preferred fit without 
adiabatic contraction, the best-fit T* and Xr are b° m lower 
with adiabatic contraction. 

6.2. Fits with constraints 

SPS models. — Table 2 gives the observed or expected B-R 
colors and the predicted using the formula log(T*) = 
-0.820 + 0.85 l(B-R) (Table 1, Bell & de Jong 2001). We 
find that the galaxy colors and the SPS-predicted values 
are significantly lower than the rotation curve-derived maxi- 
mal disk values for five out of six galaxies, the exception be- 
ing NGC 5585, although the latter is dark matter dominated 
at larger radii. In terms of velocities (Vdisk/Vtot) 2 2 ~ 0.6 for 
all six galaxies. As we have not included internal extinction 
effects, the intrinsic colors will be bluer and the resulting 



and disk velocity fractions will be lower. The maximal disk 
for NGC 5585 is inconsistent with a ~ 1 halos and cannot 
be reconciled with any halo that is adiabatically contracted. 

Sub-maximal disks. — NGC 3198 has a measured veloc- 
ity dispersion, (Vf)^ 2 = 40 ± 0.7 (Bottema 1993), imply- 
ing (Vdisk/Vobs) max = 0.43 ±0.15. For NGC 2403 we apply 
the sub-maximal disk constraint (Eq. 10). The associated 
are given in Table 2 and are consistent with the expected 
from SPS models. These constraints rule out zero and maxi- 
mal disks but do not break the disk-halo degeneracy for NGC 
2403 and NGC 3198, even when the halos are adiabatically 
contracted. 

The C-2-V200 comparison. — In Figure 10 we compare the 
fitted c-2 and V200 for a = 1 and a = halo fits with the c_2 
- V200 model of Bullock et al. (2001). This shows that as 
increases, the fitted c_2 decreases and V200 increases. For a 
given T* , fits with a = have larger c_2 and lower V200 than 
fits with a = 1 , but the differences are not significant enough to 
distinguish between a = and a = 1 . Fits with adiabatic con- 
traction follow the same path in the c_2 - V200 space as those 
without, such that for a given the c_2 of the adiabatic con- 
traction halos are lower. (Note that the c_2 and V200 values 
of the contracted halos are those of the initial pre-contracted 
halo.) For NGC 2403, NGC 3198, and UGC 2259 the Tf val- 
ues from SPS models fall within ler of the mean c_2 for a = 1 
halos, while for IC 2574, NGC 3109, and NGC 5585 the 
from SPS models require very low c_2. 

Lensing-TF. — This constraint (Eq. 14) corresponds to V200 
= 76 +21 and 87 +25 kms" 1 for NGC 2403 and NGC 3198, re- 
spectively. Note that there is no lower limit on V^oo- From 
Fig. 7, we see that for NGC 2403 the lowest V 2 m = 75 km s" 1 , 
while for NGC 3198 the lowest V 20 o = 100 kms" 1 . These val- 
ues occur when c_2 is highest, which corresponds to a = 
and T* = 0. As a and T* increase, c_2 decreases and V200 in- 
creases, so that the highest (and thus excluded) values of V200 
occur for cuspy halos with close to maximal disks. Thus, this 
constraint favors low a and fits. 

All constraints. — Our full set of constraints consist of 
from SPS models, 6 < c_2 < 30, and V200 < V max . The rota- 
tion curves for NGC 3109 and IC 2574 are still rising at the 
last data point, so we adopt V max = 100, consistent with other 
dwarf galaxies. In Table 3 we give the ranges in a that are 
within 3.5 of the best-fit x 2 - While we cannot place a sta- 
tistical confidence on these ranges, they seem to encompass 
all of our acceptable fits. Clearly, even with the full set of 
constraints significant degeneracies remain, often with both 
a ~ and a ~ 1 halos providing reasonable fits. The differ- 
ences between fits with and without adiabatic contraction can 
be substantial. Given that other evolutionary processes exist 
and are not clearly understood, there remains significant un- 
certainty in determining a. 

7. OTHER UNCERTAINTIES IN MASS MODELS 

We now discuss the sensitivity of our results to several pa- 
rameters that have so far been kept fixed. These include the 
minimum rotation curve error values, dv m i n ; the distance, D; 
the disk thickness, q^; and the halo flattening, q^. Results for 
the tests we perform below are shown in Fig. 1 1 for a = and 
a = 1 halos without adiabatic contraction. The results are sim- 
ilar with adiabatic contraction, but the range in is lower. 

Minimum rotation curve error values. — We use minimum 
error values for each velocity measurement at a given radius 
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TABLE 2 

Mass-to-light ratios from constraints 



Galaxy T« max (B-R) IWJX (B-R) obs T* 8p8 T* 8ub (VdUk/Vw)^ (W^tot)* 

(1) (2) (3) (4) (5) (6) (7) (8) 



NGC3109 1.5 1.17 0.9±0.1" 0.88™|* - 0.33 0.23 

IC2574 1.5 1.17 0.90* 0.88 - 0.41 0.24 

UGC2259 8.0 2.02 0.9±0.1 a 0-88^^ - 0.38 0.33 

NGC5585 0.7 0.78 0.86 c 0.82 - 0.47 0.35 

NGC2403 3.1 1.54 1.17±0.15 rf 1.5^;| U^i 0.50 0.42 

NGC3198 4.2 1.70 1.17±0.15 rf 1.5™ : J 1.1^ 0.59 0.56 



NOTE. — Col. (1): Galaxy name. Col. (2): Maximum disk YJ| from fits including a dark halo; for disk only models the are up to ~ 25% higher. Col. 
(3): Predicted B-R color for the in col. (2) based on SPS models from Bell & de Jong (2001). Col. (4): Observed or expected B-R color. Col. (5) 
predicted for the B—R color given in col. (4) using Bell & de Jong (2001); Col. (6) based on the sub-maximal disk constraints from Bottema (1993); 
Col. (7) ratio of disk velocity to total observed rotation velocity at 2.2 disk scale lengths for the best fit with constraints. Col. (8) ratio of maximum disk to 
maximum observed rotation velocities for the from Col. (7). 

"Expected B—R color for dwarf disk galaxies (e.g. van den Bosch & Swaters [2001] find (B—R) = 0.87 ± 0.09 for six late-type dwarf galaxies). 
b Martimbeau et al. (1994) corrected for Galactic reddening. 
c C6te et al. (1991) corrected for Galactic reddening. 

d Mean B-R color of 40 late-type HSB disk galaxies from MacArfhur et al. (2003). 




1.6 1.8 2 2.2 1.6 1.8 2 2.2 1.6 1.8 2 2.2 

log V 200 

FIG. 10. — Comparison between our fitted c_2 and V200 against the model from Bullock et al. (2001). The mean is given by the solid line, while the 1 and 2 a 
deviations are given by the long-dashed and dotted lines, respectively. Fits with a = are in black, while fits with a = 1 are in red. Fits with adiabatic contraction 
(AC) are given by filled circles and solid lines, while fits without AC (NAC) are given by open circles and dotted lines. For NGC 3109, IC 2574, and NGC 5585 
the symbols are at intervals of 0.5, while for UGC 2259, NGC 2403, and NGC 3198 the interval is 1.0. 
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TABLE 3 
Range and best-fit values of a. 



Galaxy 


OtNAC 


"AC 


a NAC+ 




(1) 


(2) 


(3) 


(4) 


(5) 



NGC3109.. 

IC 2574.. 
UGC 2259.. 
NGC 5585.. 
NGC 2403.. 
NGC 3198.. 



0-0 oio 

Q 14 
nn 



o-o ii 
o-o H 

o";i 

o.o H 



0.0° 
0.0 
0.0 
0.5 
1.2 



0.8 



03 

o.o 



0-0 : 



1.0 



1.3 
' 0.0 

0.0 ' I! 



NOTE. — Col. (1): galaxy name. Col. (2): best fit and upper and 
lower limits based on a A\l = 3.5 for fits without adiabatic contraction 
(AC) or constraints (+). Col. (3): Same as col. (2), but for fits with AC. 
Col. (4): Same as col. (4), but for fits with constraints. Col. (5): Same 
as col. (2), but for fits with AC and constraints. 



of 2, 3, and 4 kms" 1 . Larger error values have the effect of 
lowering the x 2 an d make it harder to distinguish between 
models. 

Distance. — We show the effect of a 30% distance error, in 
order to account for uncertainties in distance estimation such 
as peculiar velocities. The mass of the stars and gas is propor- 
tional to distance, so for larger distances and a given V& S k 
is lower. 

Disk thickness. — We consider the values of qo = for a thin 
disk, the commonly-used go = 1/6, and qo = 0.25. The ef- 
fect of a thick disk is to smooth out the features of the sur- 
face brightness profile and lower the amplitude of the rotation 
curve (Fig. 12). This results in a higher maximum disk but 
does not significantly alter the relative goodness of fit between 
a = and a = 1 halos. 

Halo flattening. — We compute fits for oblate halos with 
q = 0.25 and q = 0.5 and a prolate halo with q = 1.5. The 
effect of halo flattening on x 2 is very small, and is only notice- 
able for low Y* fits because the halo parameters adjust; as q 
decreases, c_2 decreases and V200 stays roughly constant. Be- 
cause we define R200 to be independent of q, then for a given 
f-2, V200 increases as q decreases. Therefore, if V200 is con- 
stant, c-2 will decrease as q decreases. Thus, although halo 
flattening does not alter the Xr> it contributes an additional 
~ 20% uncertainty to c_2. 

8. SUMMARY 

There has been much debate recently over the shape of 
galaxy density profiles, especially regarding the center of 
dwarf and LSB galaxies that are believed to be dark matter 
dominated. However, a large number of systematic effects 
such as slit position error, poorly sampled velocity fields, and 
non-circular motions thwart the straightforward interpretation 
of observed rotation curves as circular velocity profiles. 

Even if the circular velocity profile is measured perfectly, 
the determination of a is complicated by the degeneracies in- 
herent to the mass modeling exercise. The most cited degen- 
eracy is that of the unknown value of the stellar mass-to-light 
ratio, Yd; strong covariances with the halo concentration and 
density profile slope prevent a definitive determination of Y<j. 
Even if Yd were known, degeneracies that exist between the 
halo parameters might also prevent a unique determination of 
a. 

Independent constraints may help in breaking these de- 
generacies. We have considered such constraints with six 
disk galaxies that have Ha and Hi rotation curves and R- 



band imaging. The Ha rotation curves are derived from 
two-dimensional Fabry-Perot velocity fields (Blais-Ouellette 
2000; Corradi et al. 1991) and as such are not affected by 
most of the systematic caveats raised by Swaters et al. (2003a) 
on the context of long-slit spectroscopy and low-resolution 
radio-synthesis mapping, such as slit position error and beam 
smearing. However, we cannot exclude the possibility of non- 
circular motion effects in these rotation curves. 

The advantages of two-dimensional Ha velocity fields over 
long-slit spectra are demonstrated with UGC 2259, also stud- 
ied by Swaters et al. (2003a). Sampling more of the ve- 
locity field yields a scatter in the rotation curve data of 
Blais-Ouellette et al. (2004) that is significantly smaller than 
that of the long-slit spectrum of the same object by Swaters 
et al. (2003a). These authors find a = 0.86 ± 0.18 for a min- 
imum disk, and their plot of x 2 versus a is mostly flat for 
< a < 1 for all Yd. By contrast, we find the best fit a = 
for all Y*, and the x 2 increases with a. It should be noted 
that a = 1 halo fits for UGC 2259 deviate most strongly in the 
central 0.5 kpc (~ 0.01/?20o) where systematic effects on the 
rotation curve are most significant. 

For galaxies with appreciable baryonic components, the 
formation of the disk may have altered the initial dark mat- 
ter density profile. To first order the halo contracts, but other 
mechanisms such as stellar feedback and stellar bars may re- 
sult in less concentrated halos. To encompass the full range, 
we run fits with and without adiabatic contraction. Adiabatic 
contraction has the effect of turning cores into cusps, even for 
relatively low mass disks. The effect of adiabatic contraction 
on the circular velocity and density slope increases for larger 
Y*, lower c_2, and lower a. Obviously, maximal disks are 
inconsistent with adiabatic contraction. 

Applying the SPS model of Bell & de Jong (2001) to 
the expected or observed B-R colors for these galaxies im- 
plies that all galaxies are sub-maximal at 2.2Rd- This is in 
agreement with other independent techniques that suggest that 
HSB disk galaxies have, on average, sub-maximal disks with 
V d ™J. x ~ 0.6V o ^ x (e.g., Bottema 1993, 1997; Courteau & Rix 
1999; Trott & Webster 2002; Courteau et al. 2005). 

In the model from N-body simulations of Bullock 
et al. (2001) for halos with 50 < V 20 o ^ 160 kms -1 , the 2-a 
range in concentrations is 6 < c_2 < 30. By defining the con- 
centration parameter as c_2 = ?"2ooA-2, where r_2 is the radius 
where the density slope of the halo is -2, these constraints 
can be applied to halos of arbitrary a. All fits have c_2 < 30, 
but often fits with a > 1 have c_ 2 < 6. Applying this con- 
straint to mass models thus lowers the range of acceptable 
a. As a further constraint, we impose V200 < V max ; this con- 
straint does not significantly affect the best fits, but it does 
help to eliminate bad ones. If we impose the stronger con- 
straint V200 < Vmax/1-4 for bright galaxies (NGC 2403 and 
NGC 3198), values of a > 1 are disfavored. 

Without constraints NGC 3109 and IC 2574 strongly favor 
a ~ and low Yj, although both of these galaxies are not 
ideally suited for mass modeling studies (NGC 3109 has an 
uncertain inclination angle, and IC 2574 has a disrupted ve- 
locity field). The remaining four galaxies are consistent with a 
wide range of central density slopes < a < 1 .4 and mass-to- 
light ratios, Y*. Applying our full set of constraints reduces 
the range of acceptable a, but taking fits with and without 
adiabatic contraction as two extremes, there still remains a 
wide range of acceptable a, and only for NGC 5585 can we 
strongly distinguish between a ~ and a ~ 1 (in this case, 
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FIG. 1 1 . — Effect of minimum rotation curve errors, distance, disk thickness, and halo flattening for a =0 and 1 halos for all six galaxies. Minimum RC errors: 
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FIG. 11. — Continued. 
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FIG. 12. — Effect of disk thickness on disk circular velocity. The rotation curves were calculated from the observed surface brightness profile, with a disk 
thickness q& = (dashed line), 0.17 (dotted line), and 0.25 (solid line). For comparison we also show the rotation curves of exponential disks (long-dashed line) 
truncated at the last surface brightness data point and continued to infinity. The truncation radius is indicated by an arrow. The bottom panel shows the residuals 
with respect to the thin disk. 



the best-fit a = 0). Our best-fit models with constraints fa- 
vor sub-maximal disk models, with Vdisk/Vtot ~ 0-6 at 2.2 disk 
scale lengths for all six galaxies. 

Accurately determining the error bars (both statistical and 
systematic) on the observed rotation curve(s) is crucial to 
breaking the degeneracies. Doubling the minimum rotation 
curve error values from, say, 2 to just 4kms _1 reduces the 
differences in \r between models with different a or T* to 
statistically insignificant levels. 

Changing the distance or disk thickness of the galaxy alters 
the best-fitting , although the relative difference in good- 
ness of fit between a = and a = 1 halos is not significant. 
The effect of halo flattening is to decrease its concentration, 
but the effect on \ 2 is practically unchanged. 

Thus, given the above uncertainties, we conclude that ro- 
tation curve mass modeling of disk galaxies fails to provide 
tight constraints on the central density slope of dark mat- 
ter halos. 7 Constraints on central density slopes are possibly 
strongest in low-mass galaxies, especially LSB galaxies, pro- 
vided that there are no systematic errors in the rotation curves. 
Unfortunately, at present, the predictions of numerical simu- 

7 Similar limitations for the mass modeling of dwarf and LSB disk galax- 
ies are addressed in Swaters et al. (2003), and for elliptical galaxies in Cote 
et al. (2003). 



lations for these galaxy types are the weakest. 

However, the prospects for determining the relative 
amounts of dark and visible matter in disk galaxies (e.g., be- 
yond 27?d, where the rotation curve becomes flatter) look more 
promising provided that near-IR imaging and SPS models or 
velocity dispersion measurements are available to constrain 
T d . 
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